ASTROPHYSICAL JOURNAL, to be submitted 
Preprint typeset using LAT^X style emulatcapj 



STRANGE CEPHEIDS AND RR LYRAE 

J. Robert Buchler 1 , Zoltan Kollath 2 
ASTROPHYSICAL JOURNAL, to be submitted 



O 
O 

(N 
C 

a 

(N 
(N 



ABSTRACT 

Strange modes can occur in radiative classical Cepheids and RR Lyrae models. These are vibrational 
modes that are trapped near the surface as a result of a 'potential barrier' caused by the sharp hydrogen 
partial ionization region. Typically the modal number of the strange mode falls between the 7th and 12th 
overtone, depending on the astrophysical parameters of the equilibrium stellar models (L, M, T c g, X, Z). 
Interestingly these modes can be linearly unstable outside the usual instability strip, in which case they 
should be observable as new kinds of variable stars, 'strange Cepheids' or 'strange RR Lyrae' stars. The 
present paper reexamines the linear stability properties of the strange modes by taking into account the 
effects of an isothermal atmosphere, and of turbulent convection. It is found that the linear vibrational 
instability of the strange modes is resistant to both of these effects. Nonlinear hydrodynamic calculations 
indicate that the pulsation amplitude of these modes is likely to saturate at the millimagnitude level. 
These modes should therefore be detectable albeit not without effort. 
Subject headings: oscillations of stars - Cepheids - s-Cephcids 
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1. INTRODUCTION 

Strange modes were discovered and named by Saio, 
Wheeler & Cox (198) in highly nonadiabatic stars. (For a 
general reference cf. Gautschy & Saio 1996). In a recent 
paper (Buchler, Yecko & Kollath 1997; hereafter Paper I) 
it was demonstrated that strange, or surface modes not 
only exist also in classical Cepheids, thrive there. Paper I 
showed that the strange modes are vibrational modes that 
are predominantly surface modes that exist quite naturally 
even in the adiabatic limit. 

This was exhibited with the help of simple transforma- 
tions of the radial coordinate to an acoustic depth z, de- 
fined through dz — dr/c s and of the radial displacement 
vector ip — ^/pc s r8r, where c s is the sound speed, with 
which the adiabatic radial pulsation problem for radiative 
models can be reduced to a time-independent Schrodinger 
equation 
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incorporates thus the spatial variations of density p, sound 
speed c s and spherical effects, as well as gravity. This 
problem is very similar to that encountered in the study 
of wind instruments, such as the trumpet, where Eq. 1 is 
known as the Bernoulli- Webster equation. 

In Cepheids and RR Lyrae the hydrogen partial ioniza- 
tion, in which c s varies extremely rapidly, produces a very 
sharp and enormously high potential barrier, that splits 
the pulsating envelope into two almost disjoint parts (cf. 
Paper I). Using the scaled eigenvectors one readily sees 
that the 'normal' vibrational modes have a large amplitude 



in the interior region. It is however also possible for some 
modes to be predominantly trapped in the outer region - 
these surface modes are the 'strange' modes. 

When, as is customary, the modes are classified with in- 
creasing number of nodes (and thus increasing frequency) 
the strange modes appear when a quarter wavelength fits 
approximately into the surface region. For a typical 5M@ 
Cepheid on the left of the instability strip one finds that 
the strange modes have n = 10 - 12, and somewhat smaller 
n on the cold side of the instability strip because there is 
more matter outside the partial ionization region of H. 

When nonadiabatic effects, which are weak in Cepheids, 
are properly taken into account, the modal frequencies ac- 
quire of course an imaginary component. Most of the stel- 
lar envelope is damping, but the partial ionization zones 
cause strong pulsational driving. Since the strange mode 
has at best a very small amplitude in the He partial ion- 
ization region it thus experiences no driving or damping 
there. All its driving must occur in the H partial ioniza- 
tion region. It is not possible to extend the simplified de- 
scription of driving found in textbooks (e.g., Cox & Giuli 
1965) to modes higher than the fundamental or the first 
overtone. For the higher modes the driving or damping are 
a delicate function of the relative phase of the pressure and 
specific volume eigenfunctions which exhibit an increasing 
number of oscillations (related to the nodes in the adia- 
batic limit) (e.g., Glasner & Buchler 1993, Buchler, Yecko, 
Kollath 1997) as the work-integral shows 
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Fig. 1. - Cepheid model, M=8, L=10000, T cff =6700K model; 
Top: work-integrand w vs. Logio (pressure) [cgs] for modes 
(fundamental), 1 (first overtone) and 9 (strange mode); mid- 
dle: radial displacement eigenvectors for modes 0, 1, 8 and 9; 
bottom: transfored radial displacement eigenvector (^) as a 
function of acoustic depth z (for clarity mode 1 has been scaled 
up by a factor of 2 and modes 8 -10 by 12). 

In Fig. 1 (middle) wc exhibit the linear displacement 
eigenvectors (modulus) as a function of Log p of the fun- 
damental, first overtone and strange 9th overtone for a 
Cepheid model with M=8M©, L=1OOOOZ , T cff =6700K, 
X=0.70, Z=0.02. Clearly that the strange mode is very 
strongly restricted to the surface region. 

This even more apparent in the bottom of Fig. 1 which 
displays the scaled displacement eigenvectors tp ( c /- Eq. 1) 



as a function of acoustical depth z for the fundamental, 
first overtone, and for overtones 8 through 10. As pointed 
out in Paper I, this representation of the modes as ip{ z ) 
brings out the vibrating string nature of the modes, and 
shows that the strange mode (9th overtone) is mainly a 
surface mode in contrast to its neighbors. 

In the top Fig. 1 which displays the work-integrand as a 
function of Log p, the integrand is normalized so that the 
area under the curve represents the relative growth-rate 
-q = 2nkPk- The solid line represents the strange mode 
(9th overtone). For comparison we also show the work- 
integrands (multiplied by 5 for improved visibility) for the 
fundamental (dotted) and first overtone (dashed). The 
latter modes are very stable for this model which is on the 
left of the instability strip {cf. Fig. 1). The strange mode 
is restricted to the surface regions and its driving comes 
entirely from the H partial ionization region around Log p 
~3. 

Some simple a posteriori explanation of the odd stability 
properties of the strange modes can be given. The strange 
modes penetrate much less deeply than the normal modes. 
For the strange modes Wu is thus (a) totally dominated 
by the outer regions, in particular the hydrogen partial 
ionization region, and (b) Wk is much smaller. However, 
its the moment of inertia Ik is also very small because the 
mode is localized in the surface. Since the smaller work- 
integrals Wk are normalized by the smaller Ik , the strange 
modes still have relative growth-rates that are comparable 
in magnitude to those of the normal modes. However, be- 
cause only hydrogen driving contributes, the growth-rates 
of the strange modes are egregious as was already shown 
in Paper I. In that paper strange modes were found to 
be vibrationally unstable to the left (blue) of the usual 
fundamental and first overtone instability strips. 

Furthermore nonlinear hydrodynamical modelling indi- 
cated that the strange modes can give limit cycles, albeit 
with rather small pulsation amplitudes, typically in milli- 
magnitudes. This has opened the exciting possibility of ob- 
serving 'strange Cepheids', and of using such additional in- 
formation to put constraints on theoretical Cepheid mod- 
elling. 

Two questions were left unanswered in Paper I, how- 
ever. The first concerned the surface boundary condition, 
namely the disregard of the surface impedance, due to the 
reaction of the isothermal atmosphere and of running wave 
or acoustic damping. The second resulted from the neglect 
of turbulence and convection which can also affect driving 
and damping. Here, in §2 we reexamine and improve on 
the outer boundary condition. In §3 and §4 we examine 
turbulent convective Cepheid and RR Lyrae models. In 
§5 we present nonlinear pulsations of turbulent convective 
models and we conclude in §6. 

2. ISOTHERMAL ATMOSPHERE AND BOUNDARY 
CONDITIONS 

In paper I the models were integrated up to a point 
where the gas pressure was small, typically a few times the 
radiation pressure. In the linearization of the discretized 
equations we imposed a standard condition of a zero La- 
grangean pressure variation, <5p» = 0, i.e., to an open and 
perfectly reflecting boundary. This boundary condition 
disregards two effects, first, running, escaping waves, and 
second, inertial effects due to the isothermal atmosphere. 
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In acoustic parlance, these effects give rise to resistive and 
reactive surface loads (Morse & Ingard 1968). 

It is interesting to note that the change of variables from 
Sr to 'F, because of the density dependence, effectively 
changes the problem from an 'open quarter-wave length 
tube' to a 'closed half-wave length tube'. (Here we have 
taken the stellar surface to be defined as the point where 
the gas pressure equals two times the radiation pressure, 
i.e., at a very low density.) 

Another way of looking at this is to say that when lu is 
sufficiently large, the concomitant wave length gets suffi- 
ciently short compared to the the isothermal scale height 
1/a so that reflection at the surface becomes incomplete 
and acoustic energy is lost into the atmosphere. 
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Fig. 2. - Cepheid model, M=8M , L=1OOOOL . Left: growth- 
rates as a function of T e g, positive growth-rates correspond to 
linear unstable modes. The modes are labelled with the modal 
number. Right: period ratios Pk/Po as a function of T c g. The 
fundamental period Po [d] is indicated on top. The triangles 
denote the unstable overtone modes. The solid line at the bot- 
tom is the isothermal cut-off period ratio. 

In the optically thin outer region of the Cepheid model 
is essentially isothermal T — > (1/2) 1/4 T c $ . To a good 
approximation the density decreases exponentially with 
scale height a' 1 — TIT/ (fig), where g = GM/r 2 is es- 
sentially constant. This gives rise to a potential barrier 
which rapidly flattens to a constant value of V* = (c s a/2) 2 . 
In Buchler, Whiting & Kollath (2001) it is shown that 
the presence of the isothermal atmosphere can be approx- 
imately taken into account through a boundary condition 
that is imposed at the base of this atmosphere for the 
modes below the barrier (lu < lu c ). The modes above the 
barrier disappear into the continuum and are no longer 



astrophysically relevant. 

(lo + k) + (lu — k) exp(- 



S = = -ik 



■2ikA) 



(uj + k) - (lu - k) exp(-2i£;A) 



(0) 



where k = k r + i ki = yju> 2 — V*. When \ki A| is large then 
S ~ iksgnki. The approximation that is inherent in the 
expression for the boundary condition is that the poten- 
tial is a square barrier of width A. This translates into a 
boundary condition 



5p* = Tip* ( -1 + (-a/2 + S)R* Sr/R 



(7) 



With this new boundary condition our eigenvalue problem 
becomes more complicated and delicate. We compute the 
eigenvalues with a modified Castor method in which the 
Sp*, of Eq. replaces the usual 5p* — boundary condi- 
tion. The presence of strange modes and their frequent 
proximity to the regular modes makes the Castor search 
difficult. We therefore first compute all the nonadiabatic 
eigenvalues for the 5p* = case with a general eigenvalue 
solver (e.g., Glasner & Buchler 1993) and then switch on 
the new boundary condition gradually. 

P„=3.00 3.99 5.46 7.79 11.7 19.5 
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Fig. 3. 
legend. 



Cepheid model, M=6M Q , L=3OOOL ; cf. Fig. 2 for 



3. TURBULENT CONVECTIVE CEPHEID MODELS 

Turbulent convection changes the structure of the outer 
envelope, as well as the pulsational driving. To study 
its effects on the strange modes and to ascertain that 
they survive convection, we have used the Florida code 
(Yecko, Kollath & Buchler 1998) with the a's defined as 
in Feuchtinger, Buchler & Kollath (2000) with values that 
are very similar to those used in our study of first overtone 
Cepheids, viz. ctd= 2.177, a c = 0.4, a s — 0.433, a v = 0.12, 
a t = 0.001, a r = 0.4, a p = and a x = 1.5. We use the OPAL 



opacities of Iglesias 
Ferguson (1994). 



Rogers (1996) and of Alexander 
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In Figs. 2, 3 and 4 we show the results for 3 sequences of 
Cepheid models of fixed mass and luminosity over a broad 
range of T c g. The left subfigures show the growth-rates 
as a function of the T e g of the equilibrium models. The 
growth-rates for the unstable modes are positive. Some 
modes are very stable and lie below the graph. 

As a function of T c g the modal number (ordering with 
period) of the strange mode changes since it switches its 
properties with the regular modes through avoided cross- 
ings (cf. Paper I). It is therefore ambiguous how we should 
label the modes along a sequence, but this labelling is 
not important for us since we are interested in finding out 
whether there is an unstable strange mode and what its 
period is, both of which are unambiguous. 

The instability strip of the 8M© in Fig. 2 extends from 
the fundamental blue edge at T e g= 5700 K to the funda- 
mental red edge at T e $ = 4600 K. An unstable strange mode 
exists from 7700 K through 6100 K which is the 11th over- 
tone at the lower T e g, but which gradually morphs into 
the 8th overtone by the time it gets to 8000 K. A corre- 
spondingly erratic behavior shows up in the period ratios. 
It is not a numerical artifact, but arises from avoided level 
crossings (cf. Paper I) that are also familiar from the quna- 
tum mechanical behavior of an atom in a B field and from 
the nonradial p + g spectrum in unevolved stars. 



6950 K the 12th overtone is unstable, 
5100 K it is the 9th overtone. 



and from 5200 K to 



P„= 1.43 1.90 2.61 3.76 5.93 9.: 
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Fig. 4. 
legend. 



Cepheid model, M=4M Q , L=86OL ; cf. Fig. 2 for 



The instability strip of the 6M Q model in Fig. 3 ex- 
tends from the first overtone blue edge at T c s— 6200 K to 
the fundamental red edge at T e g= 5150 K. From 7000 K to 
6450 K the 10th overtone is unstable. Around 5000 K the 
8th overtone just misses being unstable. 

The instability strip of the 4M Q model Fig. 4 extends 
from the first overtone blue edge at T c g= 6550 K to the fun- 
damental red edge at T c g= 5500 K, and it includes a small 
region with second overtone instability. From 7700 K to 



4. TURBULENT CONVECTIVE RR LYRAE MODELS 

We consider an RR Lyrae model with M=O.65M , 
L=45L with X=0.70, Z=0.02. Fig. 5 shows that the in- 
stability strip of the model extends from the first overtone 
blue edge at T c g= 7250 K to the fundamental red edge at 
T cS = 6000 K. From 8300 K to 7200 K the 10th overtone is 
unstable, and from 5900 K to 5800 K the 8th overtone is 
unstable. 

Again, these results show that an excited strange mode 
can exist on either side of the classical instability strip. 

5. NONLINEAR BEHAVIOR 

The full amplitude pulsations have been computed with 
the Florida hydrodynamics code (Kollath, Beaulieu, Buch- 
ler & Yecko, 1998) with the same mesh and parameters as 
used for the linear properties. The pseudo-viscosity was 
taken to be zero in our calculations. 

P„= 0.27 0.35 0.47 0.64 0.95 
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Fig. 5. 
legend. 



RR Lyrae model, M=O.65M , L=45L ; cf. Fig. 2 for 



In Fig. 6 we exhibit the light and radial velocity curves 
for three selected Cepheid models and on RR Lyrae model. 
The lightcurves are measured in millimags, and the radial 
velocities in units of 100 m/s. Note that five of the mod- 
els (1, 3-6) are blue strange variables, i.e., they lie to 
left of the blue edge, and two (2 and 7) are red strange 
variables, i.e., they lie to right of the blue edge. Interest- 
ingly they all have about the same pulsation amplitudes. 
Because the oscillation amplitude is so small and it is so 
confined to the surface region, there is not much matter 
below the photosphere that can modulate the heat flow. 
The light-curve amplitudes are therefore also small, in the 
millimagnitudc regime, but they are observable. 

The temporal structure of the photospheric velocities 
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appears strange at first. In Fig. 7 we therefore show the 
behavior of the velocities throughout most of the model, 
in fact Lagrangean zones 65 through 145 of a 150 mass- 
shell model. The interior zones are essentially quiescent 
because of the high modal number of this mode. On the 
other hand the amplitude of the very outer zones increases 
considerably (but is not relevant because it contains so lit- 
tle mass and is considerably above the photosphere. The 
photosphere itself, which is shown as a single dashed line, is 
relatively quiescent. For reference we also show the hydro- 
gen partial ionization region with three adjacent dashed 
lines. The helium ionization lies below in zones that are 
not shown because the velocities are small. The photo- 
sphere is seen to lie in a relatively quiescent region and is 
subjected to inward and outwards moving ripples that ac- 
count for its sometimes strange looking behavior in Fig. 6. 



Lightcurve Radial Velocity 




Fig. 6. - Left: Light curves (in millimag): Right: radial ve- 
locity curves (in lOOm/s). From top to bottom, Cepheids: 
M=4M S [A], T cff =7300 and 5150 K, M=6M [C], T cff =6700K, 
M=8 [E], T cff =7100 and 6700 K, RR Lyrae: M=0.65M S [R], 
T cff =7800 and 5900 K. 

There are two shortcomings to these calculations. First, 
the (linear) special boundary condition at the outer surface 
(Eq. [?]) has not been incorporated in the nonlinear code, 
but we do not expect that to cause a major distortion, for 
the pressure is so small at the boundary. Second, what 
is perhaps more severe, there is no provision for acous- 
tic leakage in the nonlinear hydrodynamics code for waves 



with frequencies above the isothermal threshold u> c . 

The presence of higher frequency reflected waves is per- 
haps the reason for the pronounced harmonic wiggles on 
the light curves and especially in the radial velocity curves. 
Indeed, because spherical effects arc negligible, the higher 
order strange modes are all in resonance. However, if they 
are excited along with the lowest strange mode, then the 
even harmonics should be absent, but they are not. 

6. DISCUSSION AND CONCLUSION 

The improved boundary conditions have a small effect 
on the behavior of the deep modes, mostly because the 
pressure variation at the stellar surface 5p* is always small. 
However, they do affect the higher overtones as they ap- 
proach the cutoff frequency by shifting their periods, some- 
times causing inversions in the modal ordering and affect- 
ing the stability of the modes. 

Turbulent convection does not kill the linear instability 
of the strange modes, as one might have feared. On the 
contrary, they are as unstable as with the radiative mod- 
els. One even finds unstable strange modes to the right of 
the fundamental red edge. 

The growth-rates of the strange modes on the hot side of 
the instability strip are not very sensitive to turbulent con- 
vection. On the other hand convection is very strong in the 
cooler models, and the values of the stability coefficients 
of the strange modes on the cool side are consequently less 
certain than they are on the hot side. 




PHASE 



Fig. 7. - Radial velocity profile for the Cepheid model 
M—6Mq, T e flF=6700K; the individual Lagrangean zones (zones 
65 through 145 for a 150 zone model) are shifted vertically for 
better visualization; the dashed line corresponds to the pho- 
tospheric velocity, the three adjacent dashed lines denote the 
hydrogen partial ionization region 
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So, do strange Cepheids really exist in nature? We find 
the occurrence of self-excited strange modes to be quite 
robust, i.e., largely independent of any modelling param- 
eters. For an individual model the change in modelling 
parameters may have a large effect, but when one consid- 
ers a whole sequence of models on finds that modelling 
parameters (e.g., the fineness of the grid) merely cause a 
small shift (e.g., in T e g). Especially in the high luminosity 
Cepheid models, the strange mode has a period well above 
the isothermal cutoff period as the right hand side of Fig. 2 
shows. For the low mass Cepheids and the RR Lyrae 
some of the unstable strange modes occur very close to 
the isothermal cutoff frequency uj c . This cutoff frequency 
however is somewhat uncertain because it depends on the 
structure of the dilute atmosphere which itself is sensitive 
to the treatment of radiative transport (treated here in the 
customary diffusion approximation). 

We find furthermore that the full amplitude (nonlin- 
ear) pulsations are quite robust in all cases. However, 
the values of the pulsation amplitudes are somewhat un- 



certain because they depend on modelling parameters as 
they did for the classical Cepheid pulsations. In partic- 
ular the turbulent eddy viscosity has a large effect. (We 
recall that our convective parameters a's were calibrated 
to give good agreement with the observations of classical 
Cepheids, Feuchtinger et al. 2000). 

The detection of this type of pulsation would be very im- 
portant because it would add very strong additional con- 
straints on the Cepheid and RR Lyrae modelling. The 
search for millimagnitude variability on either side of the 
classical instability strip with periods approximately 1/4 
to l/5th the fundamental period (cf. Figs. 1 -3, 5) would 
be required. Hopefully this will become possible with 
COROT. 
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